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Abstract 

Using data obtained from first-principles calculations, we show that the position of the mor- 
photropic phase boundary (MPB) and transition temperature at MPB in ferroelectric perovskite 
solutions can be predicted with quantitative accuracy from the properties of the constituent cations. 
We hnd that the mole fraction of PbTi0 3 at MPB in Pb(B / B // )0 3 -PbTi0 3 , BiB0 3 -PbTi0 3 and 
Bi(B'B")0 3 -PbTi0 3 exhibits a linear dependence on the ionic size (tolerance factor) and the ionic 
displacements of the B-cations as found by density functional theory calculations. This dependence 
is due to competition between the local repulsion and A-cation displacement alignment interactions. 
Inclusion of first-principles displacement data also allows accurate prediction of transiton temper- 
atures at the MPB. The obtained structure-property correlations are used to predict morphotropic 
phase boundaries and transition temperatures in as yet unsynthesized solid solutions. 
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I. INTRODUCTION 



Perovskite solid solutions have been used extensively in ferroelectric and piezoelectric 
applications. For the past few decades, the PbZri_2.Tia.O3 (PZT) solid solution has been the 
material of choice. The highest PZT piezoelectric response is found for compositions in the 
vicinity of the morphotropic phase boundary (MPB). In the search for better fundamental 
understanding of these complex systems and for better technological properties, many other 
Pb-based solid solutions have been investigated jl} _3,_1 U,_l Currently a large amount 

of data is available in the literature on the phase diagrams and piezoelectric properties of 
a wide variety of binary and ternary Pb-based solid solutions. Nevertheless, materials with 
better electromechanical properties and/or higher transition temperature have been obtained 
only recently due to advances in single crystal growth methods |3( and due to the alloying of 
Bi ions on the perovskite A-site 0,15,0]. The search for better materials has been difficult, 
due to the lack of guidance for predicting MPB locations and transition temperatures. This 
necessitates time consuming experimental synthesis for a range of compositions to locate the 
MPB and to determine the properties of interest. While first-principles calculations are of 



great use in elucidating the structure and interatomic interactions in known materials 
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fl__L and DFT calculations are now faster than experimental synthesis, their 

utility in predicting new ferroelectric materials is weakened by an inexact representation 
of the real material by an approximate exchange-correlation functional. This often leads 
to errors in lattice constant of one to two percent Since ferroelectric materials are 

very sensitive to pressure, even such relatively small errors in volume can lead to dramatic 
inaccuracies in computed properties. Thus, the search for new piezoelectrics is currently 
guided by empirical rules based on the competition between the preferred A-0 and B-0 
bond lengths. This is encapsulated in the tolerance factor t 

f = -^, (1) 
Rb-oV2 

where Ra-o is the sum of A and O ionic radii and Rb-o is the sum of B and O ionic radii. 
For example, the tolerance factor has been widely used in solid state chemistry to predict 
stability of the perovskite structure. (A tolerance factor lower than 0.9 or higher than 1.1 
typically makes the perovskite structure unstable due to a mismatch between preferred A-0 
and B-0 bond lengths.) Recently, in a major advance, Eitel et al. suggested an empirical 



relationship [4] between tolerance factor and the transition temperature at the MPB which 
has guided the discovery of high T c alloy BiSc03-PbTi03 (BS-PT). This empirical rule states 
that as the tolerance factor of the non-PT end member decreases, the transition temperature 
at the MPB (T C MPB ) increases. 

The discovery of high T C MPB in BSPT has motivated an intense experimental effort in 
synthesis and characterization of new Bi-based piezoelectrics It was found 



that while the t-T c MPB empirical rule captures the general trend, comparison of T C MPB of 
Bi(Mg 1/2 Ti 1/2 )0 3 -PT (BMT-PT) and Bi(Mg 1/2 Zr 1/2 )0 3 -PT (BMZ-PT) revealed significant 
deviation. Whereas the t-T c MPB relationship of Ref. 4 predicts that T C MPB of BMZ-PT would 
be 80 K higher than T C MPB of BMT-PT, experimental results show that T C MPB of BMZ-PT 



is actually lower by 125 K than T C MPB of BMT-PT 

Ref. 6 considers a qualitative correlation between the tolerance factor of the non-PT end 
member and the mole fraction of PT at the MPB (xp^ PB ). They showed that while there 
is a large scatter in xf™ for higher tolerance factor values (between 1.02 and 0.95), the 
relationship becomes more quantitative for non-PT end members with lower t values (0.95 
to 0.88). 

Although the work of Refs. 4 and 6 provide some important guidelines for the search for 
new materials, several questions remain. First, what is the underlying physical reason for 
these empirical rules? Second, can better predictors of Xp PB and T C MPB be devised based on 
the understanding of the mechanisms underlying compositional phase transitions and trends 
in Curie temperatures? In this work, we show that tolerance factor is a good predictor of 
the XpT 3 due to the relationship between the local structure driven A-B cation interactions 
and phase transitions of ferroelectric perovskite solid solutions. We also show that by taking 
into account the off-centering tendencies of the B-cations as revealed by DFT calculations, 



MPB 
PT 



a more accurate predictive relationship between properties of the constituent ions and x 
and T C MPB can be obtained to guide future synthesis efforts. 



II. METHODOLOGY 

Ab-initio density functional theory (DFT) calculations were performed for the (1- 
s)Pb(Sc 1/2 Nb 1/2 )0 3 -a;PT (PSN-PT) and (l-x)Pb(In 1/2 Nbi /2 )0 3 -xPT (PIN-PT) systems 
using 2 x 2 x 2 40-atom supercells with periodic boundary conditions at experimental vol- 
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ume The energy of the system was evaluated using a local density approximation 
exchange- correlation functional and was minimized with respect to the atomic coordi- 
nates by a quasi-Newton method with no symmetry imposed. A 2 x 2 x 2 fe-point sampling 
of the Brillouin zone was used. The calculations were done with designed non-local opti- 



mized norm conserving pseudopotentials [20, |2l|. We study five compositions with x = 0, 



0.25, 0.5, 0.75 and 1.0. For pure PSN and PIN we used a 40-atom supercell with B-cations 
in a rock-salt arrangement. Three B-cation arrangements with minimal oxygen over- and 
under-bonding [3, Q were used to study x=0.25, x=0.5 and x=0.75 compositions. 
We find that the values for cation displacement vary little among the supercells studied. 
The data for PZT, Pb(Mg 1/3 Nb2/3)0 3 -PT (PMN-PT), PbfZni /a Nb y a )q a -PT (PZN-PT) 
and Pb(Sc 2 /3Wi/ 3 )0 3 -PT systems were taken from previous studies [22! . I23I I^ . The ex- 
cellent agreement between the pair distribution functions (PDF) obtained from our relaxed 



DFT structures 
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24j and by neutron-scattering experiments 
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for PMN, 



PZT and PSW indicate that our calculations accurately portray the local structure of these 
materials. 



III. RESULTS 



A. Morphotropic Phase Boundary Locations 



The correlation between Xp™ and the tolerance factor of the end-member directly reflects 
the relationship between local interatomic interactions and compositional phase transitions. 
In PT, the tetragonal phase is strongly favored, resulting in a high c/a ratio. The tetrag- 



onality is driven by off-center distortions on the A-site [27|, |28(. Perovskites where only 
the B-site cations are ferroelectrically active (e.g., BaTi0 3 and KNb0 3 ) have rhombohedral 
ground states. DFT studies have shown that the cation distortions in PbTi0 3 are coupled; 
in the absence of Ti distortion, large Pb off-centering is unfavorable, and conversely large 
Ti distortions are unfavorable if the Pb atom is held in the center of its oxygen cage [29]. 
This indicates that a strong Pb-Ti repulsive interaction is present. The existence of such 
an interaction is not surprising, since both Pb and Ti are positively charged. Since there is 
no bonding between Pb and Ti, an overlap of Pb and Ti wavefunctions gives rise to Pauli 
repulsion only. Besides the direct Pb-Ti repulsion, bond-order conservation at the oxygen 
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atoms tends to make short Pb-0 and Ti-0 bonds to the same oxygen atom unfavorable; this 
leads to an indirect oxygen- mediated Pb-Ti repulsion j^]. Similar coupling between A-site 
lone pair and B-site distortions were also noted in a wide range of complex oxides j^] . 

For Pb-based perovskites with more than one B-cation, DFT calculations show that Pb 
atoms will try to avoid the larger B-cation and move toward the smaller one l| |3| . This 
is due to a larger wavefunction overlap and a stronger repulsion with the larger ion. This 
behavior is illustrated in Figure 1 for the 50/50 PZT solid solution, with Pb ions moving 
away from the Zr-rich faces and toward the Ti-rich faces. The stronger Pb-Zr repulsion 
results in different positions of the first peak in the Pb-Zr and Pb-Ti partial PDFs (Figure 
2a). Similarly, comparison of Pb-In and Pb-Sc partial PDFs obtained from DFT calculations 
of PSN and PIN shows that first Pb-In peak is located at larger distances than the first Pb- 
Sc peak (Figure 2b). The differences in peak positions for Zr and Ti (0.11 A) and In and Sc 
(0.05 A) are similar to the Shannon-Prewitt jyj ionic radius differences between Zr and Ti 
(0.12 A) and In and Sc (0.055 A). In a solid solution, the (100) cube faces that Pb atoms 
move toward can be classified as either high, medium or low Pb-B repulsion depending on 
the size of the B-cations located at the cube face corners. For example, in the case of the 
PZT solid solutions, DFT calculations show that 3 or 4 Zr atoms in a (100) face disfavors 
Pb displacements, leading to either a smaller magnitude of Pb distortion (thus increasing 
Pb-Zr distance and lessening the Pb-Zr repulsion) or a tilt in the distortion direction toward 
a more friendly cube face with smaller number of Zr atoms. At the MPB, the Pb-B-cation 
repulsion cost induced by collinear (100) Pb distortions is equal to the cost of disordering 
the Pb distortions (due to dipole-dipole and bonding interactions). At higher Zr content, 
Pb cations distort in a variety of low symmetry directions toward available low Zr-content 
faces. This results in a disordered structure with overall (111) polarization 26, 321. 

This mechanism implies that a solid solution of PT and a perovskite with a larger size 
B-cation such as Pb(Ybi/2Nbi/2)03 (PYN) (average B-cation size of 0.754 A versus 0.72 A 
for Zr) will have an MPB at higher PT content due to stronger Pb-B-cation repulsion and 
a large energy cost of collinear (100) displacements. Conversely, a solid solution of PT and 
a perovskite with a smaller size B-cations such as PMN (average B-cation size of 0.67 A 
versus 0.72 A for Zr) will have an MPB at lower PT content due to weaker Pb-B-cation 



repulsion. 



46( . Similarities between Bi and Pb behavior on the perovskite A-site mean that 



this mechanism will not be fundamentally changed by the presence of Bi ions on the A-site. 
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We quantify these arguments by 



AE T ^ R (x nPT ) = E T (x nPT ) - E R (x nPT ) 



(2) 



and 



AE T _ R (x nPT ) = -AE dis (x nPT ) + AEa-bOiipt) (3) 

where x n px is the fraction of the non-PT end member in the solid solution, AET-R{x nP T) is 
the energy difference between the T and R phases, Ai?di S (^nPT) is the energy cost of disorder- 
ing the cation displacements due to electrostatic and bonding interaction and AEa-b^iiPt) 
is the energy cost due to A-B repulsion incurred by collinear (100) cation displacements of 
Pb and Bi ions. The MPB is located at x^ when AE T _ R (x nPT ) is and AE dis (x^) 
= A£ , A-B(^rf P r 3 )- Since Pb displacements in the tetragonal phase of PT solid solutions 
are fairly similar (Table I), we can make a zeroth-order approximation that A£'di S (x n pT) is 
independent of composition and is constant for all systems. Although this approximation is 
very rough and ignores the much larger magnitude of Bi displacements, it turns out to be 
valid for reproducing experimental observations. 

As a first approximation, we take the A-B repulsion to be a linear function of the ionic 
size of the B-cation and the A-B interatomic distance. The dependence of the AEa-b(x) 
on ionic size is intuitive - larger ions give rise to stronger interatomic repulsions. A large 
ionic size of the B-cation usually makes B-cation displacements unfavorable (rattling cation 
model) [33I 34J], giving rise to larger indirect repulsion contributions as well. Since Pb 
displacements in the tetragonal phase and the perovskite volume are fairly constant for a 
wide variety of PT solid solutions, the Pb-B interatomic distance will be a linear function 
of the B-cation displacement. 

Analysis of Bi-B-cation distances shows that the dependence of Bi-B-cation repulsion 
on ionic size and displacement of the B-cation will not be too different from the case of 
Pb-B-cation repulsive interactions. For an ideal perovskite structure with a cube edge a 
~ 4.0 A, the A-B-cation distance is 3.46 A. While a 0.45 A Pb distortion creates a Pb-B- 
cation distance of 3.23 A (assuming no B-cation displacement), the much larger 0.8 A Bi 
distortion |35J creates a Bi-B-cation distance of 3.07 A. The difference between these two 



distances is 0.16 A which is similar to the difference between the Pb 2+ and Bi 3+ ionic radii 
(1.49 A for Pb 2+ and 1.36 A for Bi 3+ ). In other words, the smaller ionic size of Bi allows it to 
move much closer to the B-cations before a further Bi displacement becomes unfavorable due 
to a rapidly rising Bi-B-cation repulsion. Since Pb and Bi are similar chemically, we expect 
that the Bi-B-cation repulsion will exhibit a similar dependence on interatomic distance and 
B-cation ionic size. The strength of the A-B repulsion is thus affected in the same manner 
by the B-cation size and off-centering for both Pb and Bi occupation of the A-site; this 
makes a fairly accurate treatment of both cases with one equation possible. We write 

A£A-B(£nPT) = £ nPT (a + b Rz S - c D^ s ) + (1 - x nPT )(ao + b R Ti - c D Ti ) (4) 

where ao, &o and Co are constants, R^ g is the average ionic size of the B-cations in the non- 
PT end member as computed from the B-cation composition and Shannon-Prewitt ionic 
radii, and D^i and D^ g are distortion magnitudes for the Ti and B-cations of the non-PT 
end member respectively. The assumption that AEt-r is a linear function of composition is 
supported by the experimental data for the PZT solid solution j3] (Figure 3). Examination 
of the data presented in Table I shows that displacements for a particular B-cation in the 
tetragonal phase are similar among different solid solutions (Figure 4). This means that 
Db is a transferable property of the B-C>6 complex. Using the displacement magnitudes 
obtained by DFT calculations, we therefore assign Db values for Mg, In, Sc, W, Zr, Fe, Nb, 
Zn, and Ti cations (Figure 9). 

At the MPB, R and T phase energies are equal and 

AE T _ R {x) = = x nPT (a + b R% g - cqD^ s ) + (1 - x nPT )(a + b R Ti - c D Ti ) - AE dis (5) 

Solving for x nPP , setting the numerator to 1 and collecting all terms that do not depend on 
the identity of non-PT end member, we get, using new constants a\,b% and c\ 

x nPT = l/(a 1 + b 1 R^-c 1 D^). (6) 
Expanding this expression in a Taylor series to first order, we get 
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XnPT = a 2 -b 2 R™ g + c 2 D™ g , (7) 

where a 2 , b 2 and c 2 are new constants resulting from the Taylor expansion. Since the PT 
content at MPB (xf^ B ) is equal to l-a^ P p T B , we get 



a#PB = 1 _ fl2 + b 2 R^ - c 2 D^ g (8) 

Since a large average size of end member B-cations is equivalent to a lower end member 
tolerance factor, we can write 



^ T PB = Os - MnPT - C 3 D^, (9) 

where we once again use a Taylor series to convert from R p g to t n pT and collect all terms 
that do not vary with t n pT or R p g into the constant 03. 

The dependence of x^ PB on t n pT in Eqn. 9 explains the empirical relationship between 
t and xp^ B found in Ref. 6. Examination of the plot of xp^ B vs t n pT (Table II, Fig 5) 
shows that a narrow band of about 0.2 width encompasses most points, with outliers lying 
below the lower edge. (PbMn 2 /3Wi/ 3 3 is the exception to this trend, but perhaps the 
presence of the much larger Mn 2+ ions makes the tolerance factor assigned to this system 
incorrect as discussed below). For some cases the deviations are quite large, e.g. for PZN 
and PbFei/2Nbi/203 (PFN). For solid solutions with Xp PB values lying in the narrow band, 
no more than two thirds of the end member B-cations are occupied by ferroelectrically 
active Nb and Fe ions and no strongly off-centering Zn and Ti cations are present. For 
these solid solutions, the c 3 D p g term is fairly constant, and the position of the MPB shows 
a good correlation with the tolerance factor of the non-PT end member. On the other 
hand, examination of the outlier compositions reveals that the B-site is either fully occupied 
by Fe 3+ and Nb 5+ B-cations or contains Zn and Ti cations. The large Nb, Zn and Ti 
displacements found in our DFT calculations and the 0.17 A Fe 3+ displacement found by 
DFT calculations on BiFeC>3, mean that D p g values for these solutions are significantly 
higher than the typical D p g values found for the data points in the narrow band. Thus, 
the deviations from the t-xp^ B relationship are due to differences in the average B-cation 

8 



displacements. Sorting the solid solutions by the content of ferroelectrically active cations on 
the non-PbTi0 3 end-member B-site, we see that the data now falls into two distinct groups 
with the lower x^ B values for the end-members with high content of ferroelectrically active 
cations (Figure 5). 

We now quantify the above qualitative analysis. Since correlations based on Eqn. 8-9 
will hold only as long as the changes in the ionic displacements and radii are not too large 
and the Taylor expansion converges, we choose the functional form in Eqn. 5 in order to 
treat the full range of possible PT solid solutions. Using R B and D B values (Table III) for 
solid solutions containing elements for which we have DFT D B data (Table I), we fit the 
experimentally observed PT content at MPB to equation 6 to obtain a new predictor of 

™MPB . 
PT 



x mpb = 1 _ 1 /(_ .32 + 4.53i£ vg - 7.99L>* vg ). (10) 

The correlation between the xp^ B values predicted by the fit and observed experimentally 
is shown in Fig 6. All predicted PT mole fractions are within 0.15 for the end members with 
B-cations for which we have DFT displacement data. 

Examination of the correlation in Figure 6 shows that Eqn. 10 reproduces general trends 
in MPB locations for solid solutions with only Pb on the A-site as well as for Bi(B'B")03-PT 
solid solutions with both Pb and Bi ions on the A-site. This confirms that while the smaller 
ionic size of Bi decreases Bi-B-cation repulsion making the tetragonal phase more stable, 
the larger Bi displacement has the opposite effect of almost the same magnitude. 

As a check for obtained results, we evaluate xp^ B for PT end member itself. Clearly, 
mixing PT with PT will never result in an MPB. Therefore, a prediction of a positive value 
of XpT B by Eqn. 10 would indicate a failure of our model. In fact, using R Tl and D T \ 
values from Table II, we get Xp| >B =-1.47. This means that PT is tetragonal, as observed 
experimentally. 

Despite the lack of DFT displacement data, simple arguments can be used to estimate 
reasonable D B values for Co 2+ , Ni 2+ , Cd 2+ , Mn 3+ and Yb 3+ cations. These can then be 
used to further test the validity and the predictive usefulness of Eqn 10. The ionic size of 
Ni 2+ , Co 2+ , Cd 2+ and Yb 3+ (in the high-spin state typically found in perovskites) is similar 
to or larger than that of Mg 2+ cation. According to the rattling cation model, such a large 
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size will make distortions unfavorable as the large B-cations do not have room to off-center. 
Unlike the Zn 2+ cation, Ni 2+ , Co 2+ and Yb 3+ do not have empty low lying y electronic states 
that enable covalent bonding and distortion of Zn 2+ in the Zn-Og cage |23l. |27|. In the case 
of Cd 2+ (where such low lying states are present), the extremely large ionic radius (0.95 A) 
will prevent significant off-centering. We therefore estimate values for these cation to be 
equal to 0.08 A obtained for Mg 2+ by DFT calculations. The smaller Mn 3+ has more room 
to off-center and should be more ferroelectrically active with Db value higher than that for 
Mg 2+ . The lower Ti content at MPB for PFN-PT than for PbMn 1/2 Nbi /2 03-PT (PMnN- 
PT) despite the identical size of the Fe 3+ and Mn 3+ cations, implies that D^n should smaller 
than Dp e (0.17 A). Based on experimental MPB data for Fe and Mn containing systems, we 
estimate D Mn of 0.12 A. 

The x^ BB values generated for Pb(B'B")03-PT solid solutions based on the assigned 
displacements (Table III, Figure 6) are generally in good agreement with experimental val- 
ues jUll- However, large disagreement is found for Pb(Mn 2 /3Wi/3)-PT solution as well as 
a general underestimation of x^ B for BiB03-PT and BiB'B"03-PT systems and a general 
overestimation of x^ BB for Pb-based perovskites. In the case of Pb(Mn 2 / 3 Wi/3)-PT it is pos- 
sible that this discrepancy (which is also present in Figure 5) is due to the presence of Mn 2+ 
ion in the system. While the difference between Pb(Mn 2 / 3 Wi/3)03 and Pb(Mni/ 2 Wi/ 2 )03 
Mn mole fractions is fairly small (0.16), the impact of the change of Mn 3+ into Mn 2+ is signif- 
icant because of the large difference in ionic radii of the two Mn oxidation states (0.83 A for 
Mn 2+ and 0.645 A for Mn 3+ ). The larger ionic radius of Mn 2+ gives rise to smaller tolerance 
factor and to a higher PT content at the MPB than that predicted based on Mn 3+ ionic 
size. In fact, 0.55 MPB PT mole fraction predicted for the PbMni/ 2 Wi/ 2 03-PT solution by 
Eqn. 10 is quite close to the experimental value of 0.53. 

The underestimation of x^ B for BiBC>3-PT and BiB'B"03-PT systems suggests that Bi 
substitution on the A-site gives rise to a change in dependence of the A^a-b on the average 
ionic size and displacement of the B-cation. Although this change is not large as indicated 
by the moderate size of the deviations of predicted x^ B from experimental values, the 
effect of large Bi displacement is not fully canceled out by a smaller Bi ionic size, leading 
to stronger A-B repulsive interactions and underestimation of x^ BB and overestimation of 
for Pb-based perovskites. 

To take the difference between Pb and Bi into account, we now fit Eqn. 6 separately for 
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Pb-based solid solutions and for solutions with mixed Pb/Bi occupation on the A-site and 
obtain 

x mpb = x _ 1 /(o_ 34 + _ 7.49^ vg ). (11) 

for solid solutions with only Pb on the A-site, and 

a#£ B = 1 - l/(-0.97 + 5.71i^ vg - 7.69Z^ vg ). (12) 

for solid solutions with Pb and Bi on the A-site. This improves the correlation between the 
predicted and experimental Xp PB values (Figure 7). We therefore expect that Eqns. 11-12 
will be better predictors of values than Eqn. 10. 



B. Transition Temperatures at the MPB 

We find that incorporation of the first-principles displacement data leads to more accu- 
rate predictions of T c at the MPB. Eitel et al. have shown that despite some scatter in 
the experimental data, a lower tolerance factor of the non-PT end-member generally leads 
to a higher T C MPB . Joint experimental and DFT studies of Pb-based solid solutions have 
shown that ferroelectric transition temperatures in perovskite systems are proportional to 
the square of Pb and B-cation displacements |23L l24l| and are enhanced by coupling between 
nearest neighbor ferroelectrically active cations |38j. We therefore investigate if the scatter 
in the t-T c MPB correlation is due to the differences in B-cation displacements. 

Examination of the experimental T C MPB data (Table IV, Figure 8) shows that the data fall 
into two groups. Both groups exhibit a linear dependence of J' C MPB on t but the T C MPB val- 
ues are shifted down for BS-PT, BMZ-PT and BiMg3 /4 W 1/4 03-PT (BMW-PT), PSW-PT, 
PbCo 1/2 W 1/2 3 -PT (PCoW-PT) and PbMg 1/2 W 1/2 3 -PT (PMW-PT), relative to T C MPB 
values for the other solid solutions. This split is due to the differences in B-cation com- 
positions between the two groups. The end-member perovskites in the first group (from 
now on referred to as group 1) do not contain any strongly ferroelectric Nb, Fe, Ti or Zn 
ions on the B-site. Instead, the B-site is fully occupied by weakly ferroelectric cations Mg, 
Co, Sc and W in BS, BMW, PCoW and PSW end-members. In BMZ end-member, one 
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half of the B-site is occupied by weakly ferroelectric Mg cation and one half of the B-site is 
occupied by the Zr cation, which displays strong off-centering only in solid solutions with 
high PT concentration. On the other hand, for the end-member perovskites in the second 
group (from now on referred to as group 2) either at least half of the B-site is occupied by 
the strongly off-centering Nb, Fe, Ti or Zn ions or the absence of strongly off-centering ions 
in the non-PT end-member is balanced by the absence of the weakly ferroelectric ions in the 
case of PbZr0 3 -PT, PbHf0 3 -PT and PbSn0 3 -PT solutions 

The dependence of T C MPB values on ferroelectric activity of the B-cations is consistent 
with the correlation between Curie temperatures and the square of ionic displacements in 

l39| . Similarly, a correlation between a high T c and low t 
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24 



ferroelectric perovskites 
(which has been noted previously |4j) can be explained as follows. Small t values indicate 
that the B-site cation prefers a larger volume than the A-site cations. Such an expansion of 
the A-site has been shown to lead to increased A-site displacements for the ferroelectrically 

27 . [jO^ . Since transition temperatures have been shown 



active Pb and Bi A-site cations 



to correlate with the square of ionic displacements, a small tolerance factor will lead to a 
higher transition temperature. Especially small end-member tolerance factors are achieved 



for Bi-based perovskites, due to the sma' 



1 size of Bi 3+ cation. The larger separation of Bi 
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42| may also contribute to the high T c values 



lone pairs as compared to Pb lone-pairs 
observed for Bi-based ferroelectrics. 

Carrying out linear regression analysis of the T C MPB data for the two groups of solid 
solutions 47] we obtain 

T C MPB = 4511 - 4514t (13) 
for solid solutions in group 1, and 

T C MPB = 6634 - 6539t (14) 

for solid solutions in group 2, where T C MPB is in °C. The fit to data in group 1 is quite similar 
to the extrapolation used by by Eitel et al. to predict T C MPB in BS-PT, Biln03-PT (BI- 
PT) and BiYbCvPT (BY-PT) solid solutions. This is consistent with the weak ferroelectric 
behavior of Sc 3+ , In 3+ and Yb 3 + ions. The relatively low scatter of the data in Fig. 8 means 
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that Equations 13 and 14 can be used to accurately screen for perovskite solid solutions with 
high T C MPB values. 

IV. DISCUSSION 

The strong correlation exhibited in Figure 7 suggests that the dependence of the MPB 
position on the ionic size and displacements should hold for as yet unexamined Pb and 
Bi-based solid solutions. Below we discuss some of the solutions which should exhibit 
interesting properties at the MPB (Table V). Among Pb-based binary perovskites, the 
PbZni /2 Wi /2 3 (PZW) and PbCd 1/2 Wi /2 3 (PCdW) double perovskites and their solid 
solution with PbTiC>3 are the most promising systems for experimental investigation. Un- 
like the PbB^~ 2 B^J" 2 03 perovskites, all PbB^B^Os perovskites exhibit a rocksalt B-cation 
arrangement due to a high value of disordering energy created by large charge difference 
of +2 and +6 cations. These compounds are typically antiferroelectric; the ferroelectric 
rhombohedral phase only appears after substantial (more than 0.2) Ti substitution and is 
followed by a tetragonal phase For PZW, substitution of the highly ferroelectrically 
active Zn 2+ on the B 2+ site should have an impact similar to Ti substitution. At the very 
least, this will shift the AFE-FE transition as well as the MPB to lower x^ B value (Table 
V). It is also possible that the end member itself will be ferroelectric in the ground state 
due to total B-cation ordering and the partial occupation of the B-site by ferroelectrically 
active cations. For the PCdW-PT solid solution, Eqn.13 predicts that T 6 MPB (487 °C) will 
be higher than that of PZT (385 °C). 

For BiB'B"03-PT solid solutions, Eqn. 12 predicts that a combination of Zn on the B' 
site and Ti, Nb or W cations on the B" site will require very low PT substitution to stabilize 
the tetragonal phase. Using R^ s and D# s values in Table II, we get Xp PB values of -0.07, 
0.26 and 0.21 for BiZni/ 2 Tii/ 2 03, BiZn 2 / 3 Nbi/303 and BiZn 3 / 4 Wi/403 respectively. The 
negative value of Xp PB for BiZn^Ti^Os implies that mixing BiZn^Ti^Os with PbTiC>3 
will never result in an MPB, with tetragonal phase being preferred for all compositions. 
This prediction is consistent with our recently obtained experimental data for BZT-PT 
solid solutions [3] (Figure 9). 

One of the main goals of current research in piezoelectrics is the discovery of solid solutions 
with high T c at the MPB. For the new systems to be technologically useful, it must be 
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possible to synthesize a pure perovskite phase of the material under ambient conditions. 
While extremely low tolerance factor end-members (e.g. BiYbOs) would give rise to high 
PT content at MPB and a high T C MPB , the instability of the perovskite phase in the pure end 



member make the synthesis of such solid solutions challenging |4J, [16 , Q. BiB'B"0 3 end- 
members with high ferroelectric cation content provide an alternative route for achieving high 
MPB transition temperatures without sacrificing perovskite stability. For BiCdi/2Tii/203 
(BCdT), BiZn 1/2 Zr 1/2 3 (BZZ), BiZn 2/3 Nb 1/3 03 (BZN) and BiZng^W^Og (BZW) end- 
members, Eqn. 14 predicts T C MPB values of « 770 °C, 640 °C, 580 °C and 570 °C, respectively. 
These compare favorably with j' c MPB values of 507 °C and 615 °C for BI and BY end- members 
respectively, predicted by Eqn. 13. For BZW and BZN solid solutions, Eqn. 11 predicts a 
rather low PT content at MPB. It has been shown experimentally that for BZN-PT solid 
solution, pure perovskite phase stability extends only to 30% BZN content |44j. This will 
make the synthesis of the MPB composition difficult by conventional solid-state chemistry 
methods. For BCdT end-member, Xp PB value is 0.52. Since the BS-PT solid solution is 
stable in the perovskite phase up to 0.5 BS content, the prospects for synthesizing MPB 
compositions with very high Curie temperatures in this system are more promising. 



V. CONCLUSIONS 



Using experimental data from the literature and our own first-principles calculations, we 
have shown that the position of the MPB and T c at the MPB in PbTi03 based solid solutions 
can be predicted from the average ionic radius and the B-cation displacement of the non- 
PbTi03 end member. Inclusion of first -principles data leads to significant improvements on 
the previously observed empirical correlations, allowing quantitatively accurate predictions 
of MPB location and T C MPB . In agreement with previous work, large average B-cation radius 
(or lower tolerance factor) is found to favor high PT content at the MPB. Large values 
of B-cation displacements give rise to lower PT content at the MPB. Our model predicts 
that tetragonal phase will always be favored in BZT-PT solid solution and suggests that 
PZW-PT solid solution will exhibit a phase digram different from other Pb(B 2+ ,B 6+ )03-PT 
solutions. Based on the relationship between tolerance factor and T C MPB , we suggest that 
BiZn^Zr^Og-PT and BiCd^Ti^Os-PT solid solutions are the most promising candi- 
dates for the synthesis of high T c piezoelectrics. 
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TABLE I: Results of our DFT calculations for ferroelectric compositions of PbTi03 solid solutions. 
Ground-state cation displacements from center of oxygen cage in A. Rhombohedral, monoclinic and 
tetragonal phases are denoted by R, M and T respectively. 

Pb,Bi B' B" Ti Phase 
disp disp disp disp 



PMN-0.25PT 


0. 


.389 


0. 


.080 


0. 


162 


0. 


220 


R 


PMN-0.63PT 


0. 


.387 


0. 


099 


0. 


.174 


0. 


233 


T 


PT 


0. 


,440 










0. 


.280 


T 


PZN-0.25PT 





.461 


0, 


.258 


0. 


.171 


0. 


.272 


T 


PZN-0.63PT 


0. 


,424 


0. 


.270 


0. 


.187 


0. 


.256 


T 


PT 


0. 


.440 










0. 


.280 


T 


PSW-0.25PT 


0. 


.344 


0. 


.095 


0. 


070 


0. 


.241 


R 


PSW-0.63PT 





,350 


0. 


.087 


0. 


.126 


0. 


.208 


T 


PT 


0. 


.440 










0. 


.280 


T 


PZ-33PT 


0. 


,462 


0. 


.129 


0. 


.129 


0. 


235 


R 


PZ-50PT 


0. 


.463 


0. 


.138 


0. 


.138 


0. 


.247 


M 


PZ-66PT 


0. 


,495 


0. 


.175 


0. 


.175 


0. 


.278 


T 


PT 


0. 


.440 










0. 


.280 


T 


PSN-0.25PT 


0. 


.402 


0. 


.112 


0. 


.177 


0. 


236 


R 


PSN-0.50PT 





.381 


0. 


.112 


0. 


.155 


0. 


.243 


M 


PSN-0.75PT 


0. 


.411 


0, 


.143 


0. 


166 


0. 


.240 


T 


PT 





.440 










0, 


.280 


T 


PIN-0.25PT 


0. 


.407 


0. 


090 


0. 


.156 


0. 


292 


R 


PIN-0.50PT 





,368 


0. 


073 


0. 


.162 


0. 


.272 


T 


PIN-0.75PT 


0. 


.385 


0. 


.083 


0. 


.178 


0. 


.261 


T 


PT 





.440 










0. 


.280 


T 
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TABLE II: Position of the MPB in PbTi0 3 solid solutions. Data for BMW-PT and PSW-PT 
solutions taken from Ref. 18 and 20 respectively. Data for all other systems taken from Ref. 
6. End-members with high ferroelectric activity on the B-site are marked by an asterisk. For 
Pb(Mn,W)0 3 -PT solid solution, data is given for PbMn 2 /3,Wi/ 3 03 and PbMni/2,Wi/2)0 3 possible 
end-member stoichiometries. For PbFe2/ 3 ,Wi/ 3 3 -PT solid solution the average position of the 
MPB region that extends from 0.2 to 0.4 PT content |45j is given. 

/ T MPB 



PbFe9/-?Wi rtCh 


1.007 


0.30 


PbMno/qWl /qO'! 


1.007 


0.53 


PbFei/sNbi/sOs* 


1.001 


0.10 


PbMni /oNbi/oOs 


1.001 


0.25 


PbNii rtNbortOq 


0.994 


0.35 


PbMgi „W, /oOq 


0.992 


0.45 


PbMg 1/3 Ta2/303 


0.989 


0.32 




0.989 


0.33 


PbCd „W, /9O3 


0.986 


0.45 


PbZni /oNb<mCh* 


0.986 


0.10 


Nai /oBii /oTiOq* 


0.980 


0.14 




n QR4 


n w 

u.oo 


PbSn0 3 


0.978 


0.45 


PbSc 1/2 Nb 1/2 3 


0.977 


0.43 


PbScx^Tax^Oa 


0.977 


0.45 


PbHf0 3 


0.969 


0.52 


PbZr0 3 


0.964 


0.48 


PbSc2 /3 W 1/3 03 


0.974 


0.45 


PbIn 1/2 Nb 1/2 3 


0.964 


0.36 


Pbin 1/2 Ta 1/2 3 


0.964 


0.38 


BiFe0 3 * 


0.956 


0.28 


BiMn0 3 


0.956 


0.50 


BiNi 1/2 Ti 1/2 3 * 


0.955 


0.45 
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i „,MPB 



PbYb 1/2 Nb 1/2 3 0.949 0.50 

BiMg 1/2 Ti 1/2 3 * 0.948 0.34 

BiMg3 /4 W 1/4 03 0.934 0.50 

BiMg 1/2 Zr 1/2 3 0.922 0.60 

BiSc0 3 0.911 0.65 
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TABLE III: Ionic data and predictions for PT content at MPB using Eqns. 11-12 for Pb- and 
Bi-based end member perovskites. Shannon-Prewitt B-cation ionic radii and displacements ob- 
tained by DFT calculations are in A. Fe displacement value is taken from Ref. 36. B-cation 
displacement data marked with an asterisk are estimated based on crystal chemical arguments. 
For Pb(Mn,W)03, data is given for PbMn^W^Os and PbMn2/3W!/ 3 03 stoichiometries. 





B',B" 


t 


B 


MPB 


MPB 




size 




disp 


predict 


exp 


PMW 


0.72,0.60 


0.992 


0.08,0.10 


0.46 


0.45 


PCoW 


0.745,0.60 


0.986 


0.08*,0.10 


0.47 


0.45 


PMnW 


0.83,0.60 


0.966 


0.08*,0.10 


0.51 


0.53 


PMnN 


0.645,0.64 


1.001 


0.12*,0.17 


0.28 


0.25 


PFN 


0.645,0.64 


1.001 


0.17,0.17 


0.16 


0.10 


PSN 


0.745,0.64 


0.977 


0.11,0.17 


0.37 


0.43 


PIN 


0.80,0.64 


0.964 


0.07,0.17 


0.45 


0.36 


PYN 


0.87,0.64 


0.949 


0.08*,0.17 


0.48 


0.50 


PNN 


0.69,0.64 


0.994 


0.08*,0.17 


0.32 


0.35 


PMN 


0.72,0.64 


0.989 


0.08,0.17 


0.33 


0.33 


PZN 


0.74,0.64 


0.986 


0.25,0.17 


0.09 


0.10 


PSW 


0.745,0.60 


0.975 


0.11,0.10 


0.46 


0.45 


PFW 


0.645,0.60 


1.002 


0.17,0.10 


0.25 


0.30 


PMnW 


0.645,0.60 


1.007 


0.12*,0.10 


0.37 


0.53 


PZ 


0.72,72 


0.964 


0.13,0.13 


0.43 


0.48 


BF 


0.645,0.645 0.956 


0.17,0.17 


0.29 


0.28 


BMn 


0.645,0.645 0.956 0.12*,0.12* 


0.44 


0.50 


BS 


0.745,0.745 0.911 


0.11,0.11 


0.59 


0.65 
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B',B" t B MPB MPB 

size disp predict exp 



BNiT 0.69,0.60 0.955 
BMgT 0.72,0.60 0.948 
BMgZ 0.72,0.72 0.922 

BMgW 0.72,0.60 0.934 



0.08*,0.25 0.31 0.45 

0.08,0.25 0.35 0.34 

0.08,0.13 0.57 0.60 

0.08,0.10 0.57 0.48 
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TABLE IV: T c at the MPB in PbTi0 3 solid solutions. Data for for Bi(Mg,W)0 3 -PT and 
Pb(Sc,W)03-PT solid solutions taken from Ref. 18 and 24 respectively. Solid solutions with 
lower ferroelectric activity on the B-site are marked by an asterisk. 



Pb(Co,W)0 3 * 


0.987 


30 


Pb(Fe,W)0 3 


1.001 


52 


Pb(Mg,W)0 3 * 


0.993 


60 


Pb(Mg,TaO) 3 


0.989 


80 


Pb(Sc,W)0 3 * 


0.975 


97 


Pb(Ni,Nb)0 3 


0.994 


100 


Pb(Fe,Nb)0 3 


1.001 


130 


Pb(Mg,Nb)0 3 


0.989 


140 


Pb(Mn,Nb)0 3 


0.973 


160 


Pb(Zn,Nb)0 3 


0.986 


187 


Pb(Sc,Ta)0 3 


0.977 


190 


PbSn0 3 


0.978 


205 


Pb(Co,Nb)0 3 


0.985 


220 


Bi(Mg,W)0 3 * 


0.933 


225 


Pb(Sc,Nb)0 3 


0.977 


260 


Bi(Mg,Zr)0 3 * 


0.921 


300 


Pb(In,Nb)0 3 


0.965 


320 


PbHf0 3 


0.969 


340 


(Na,Bi)Ti0 3 


0.977 


350 


Pb(Yb,Nb)0 3 


0.951 


360 


PbZr0 3 


0.964 


385 


Bi(Ni,Ti)0 3 


0.955 


425 


Bi(Mg,Ti)0 3 


0.947 


425 


BiSc0 3 * 


0.911 


450 
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TABLE V: Predicted x^ B and T C MPB values for PbB'B"0 3 -PT and BiB'B"0 3 -PT solid solutions 
discussed in Section IV. Predicted £11*6 obtained using Eqn. 11-12, T^ PB (in °C) are obtained 

by uisng Eqn. 13-14. Ionic size and displacement data in A. 

B',B" B',B" t AvgB xf BB T C MPB 
size disp disp predict predict 



Biln0 3 
BiYb0 3 



90 




1 80 

-LOU 


0.08 


0.55 


487 


0.05 


0.67 


507 


0.05 


0.70 


615 


0.15 


0.54 


770 


0.19 


0.42 


644 


0.22 


0.26 


579 


0.21 


0.30 


572 
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5 



6 



7 



8 



FIG. 1: (Color online). Projection of the DFT relaxed structure for 4x 2 x 1 50/50 PZT supercell 
on the x — y plane. Pb atoms tend to distort aw 




Distance ( A ) Distance (A) 



FIG. 2: Pb-B-cation partial PDFs. (a) Pb-Ti (solid) and Pb-Zr (dashed) partial PDFs obtained 
from the relaxed structures of the 50/50 PZT supercells at experimental lattice constants, (b) 
Pb-Sc (solid) and Pb-In (dashed) partial PDFs obtained from the relaxed structures of PSN and 
PIN at experimental lattice constants. 
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FIG. 3: Experimental AEt-r values (circles) for the PZT solid solutions as function of Ti com- 
position 36]. The dependence of AEt-r on PT is well fit by a linear function for the whole 
compositional range. 



0.25 



§ 0.2 

g 

o 

U0.15h 

CO 

Q 

.2 0.1 

cd 
o 

! 

P3 

0.05 



Weakly FE active 
▲ Moderately FE active 
+ Strongly FE active 



: t * 



Me 



In 



A 

I 



♦ 
♦ 



Sc 



W 



Zr 



Fc 



Nb 



Zn 



t 



♦ 
♦ 



Ti 



FIG. 4: (Color online). Distortion magnitudes (in A) for B-cations in PT solid solutions. Distortion 
magnitudes are especially large for Zn and Ti. Z?b values of 0.08 A, 0.07 A, 0.11 A, 0.10 A, 0.13 A, 
0.17 A, 0.17 A, 0.25 A and 0.25 A are assigned for Mg, In, Sc, W, Zr, Fe, Nb, Zn and Ti respectively. 
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0.92 0.94 0.96 0.98 1 1.02 

t of non-PbTi0 3 end-member 

FIG. 5: (Color online). Plot of the MPB position (in mole fraction PT) versus the end member 
tolerance factor. Circles indicate non-PT end-member perovskites with no Ti-cations and less than 
full occupation of the B-site by ferroelectrically active Nb or Fe cations. These display a lower 
ferroelectric activity on the end-member B-site. Stars indicate systems where ferroelectric activity 
of the non-PT end- member B-site is higher due to presence of Ti ions and/or occupation of more 
than one half of the B-sites by ferroelectrically active Fe and Nb ions. Diamonds indicate data 
for Pb(Mn,W)03-PT solid solution with 1:1 (low tolerance factor) and 2:1 Mn/W (high tolerance 
factor) stoichiometries. 
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0.2 0.4 0.6 

Experimental x pT at MPB 



FIG. 6: Correlation between the mole fractions of PT at MPB predicted by Eqn. 10 using data 
in Table III and mole fractions of PT at MPB observed experimentally. Solid solutions for which 
we have DFT B-cation displacement data are marked by filled cirles and diamonds for Pb-based 
and Bi-based systems, respectively. MPB positions predicted for Pb- and Bi-based systems where 
B-cation displacement data is estimated are marked by open circles and diamonds respectively. 
Predicted positions for MPB for Pb(Mn,W)G"3-PT solid solutions using 1:1 (predicted Xprp B =0.55) 
and 2:1 (predicted Xp BB =0.39) Mn:W stoichiometrics are represented by stars. Taking the effect 
of B-cation off-centering into account makes quantitative prediction of MPB location possible. 
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FIG. 7: Correlation between the mole fractions of PT at MPB predicted by Eqn. 11-12 using data 
in Table III and mole fractions of PT at MPB observed experimentally. Solid solutions for which 
we have DFT B-cation displacement data are marked by filled cirles and diamonds for Pb-based 
and Bi-based systems respectively. MPB positions predicted for Pb- and Bi-based systems where 
B-cation displacement data is estimated are marked by open circles and diamonds respectively. 
Predicted positions for MPB for Pb(Mn,W)0"3-PT solid solutions using 1:1 (predicted Xpr^ B =0.51) 
and 2:1 (predicted x^ B =0.37) Mn:W stoichiometrics are represented by stars. Using separate fits 
for Pb- and Bi-based end-member perovskites improves the correlation. 
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FIG. 8: (Color online). T^ PB versus tolerance factor of the non-PT end- member using data of 
Table V. End-members with lower FE activity on the B-site are shown by circles, end-members 
with higher FE activity on the B-site are shown by squares. Low concentration of ferroelectrically 
active cations on the end-member B-site shifts the T^pb ^ q i ower values. Least-squares fit lines 
for the two cases are shown. 
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FIG. 9: Experimental data for c/a ratio in BiZn^Tii^Os-PT solid solution. The material is 
tetragonal for all compositions studied, with an anomalous strong enhancement in tetragonality. 
Compositions with less then 0.6 PT content could not be made as a single phase. 
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